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Generalizing the first law of thermodynamics, the increase in entropy density 5S{x) of a conformal 
field theory (CFT) is proportional to the increase in energy density, 6E{x), of a subsystem divided 
by a spatially dependent entanglement temperature, Te{x), a fixed parameter determined by the 
geometry of the subsystem, crossing over to thermodynamic temperature at high temperatures. In 
this manuscript we derive a generalization of the thermodynamic Clausius relation, showing that 
deformations of the CFT by marginal operators are associated with spatial temperature variations, 
6Te{x), and spatial energy correlations play the role of specific heat. Using AdS/CFT duality we 
develop a relationship between a perturbation in the local entanglement temperature of the CFT 
and the perturbation of the bulk AdS metric. In two dimensions, we demonstrate a method through 
which direct diagonalizations of the boundary quantum theory may be used to construct geometric 
perturbations of AdSs. 


Important connections have recently been forged be¬ 
tween quantum information and gravitational physics. A 
relationship between the change in entanglement entropy, 
5S^ of a subregion of a conformal field theory (CFT) and 
an energy, 5E^ associated with an excited state of the 
CFT has been developed that generalizes the first law 
of thermodynamics [1]. Subsequently it was shown that 
the resulting generalized first law (GFL) was equivalent 
to the linearized Einstein Equations for perturbations of 
the AdS bulk [2j [3] . The GEL was further elucidated by 
Klich and coworkers [4] who identified an ’’entanglement 
temperature” function by showing that the reduced den¬ 
sity matrix of a subregion in a CET could generally be 
written: 

p— ! dx5{x)f{x) ^ 2 ) 

where T{x) is the energy density operator of the CET 
and P{x) is the inverse of the entanglement temperature 
that diverges at the boundary of the subsystem. Within 
this framework, the GEL is simply derived by perturbing 
the ground state density matrix hy 5p and expanding the 
entropy, S = —trp log p to first order to yield: 

SS = J dxP{x)6T{x) (2) 

where 5T{x) = tT{5pT{x)) 

In analogy to an equation of state, it is natural to 
ask how the entropy responds to a change in an exter¬ 
nal parameter such as an applied field. To see that some 
relation other than the GEL is needed we note that one 
might easily apply a field that increases the energy while 
decreasing the entanglement entropy. A simple exam¬ 
ple, which we subsequently study, is a one dimensional 
noninteracting spinless fermion gas, corresponding to a 
c = 1 CET in the continuum. The marginal backscatter- 
ing operator locally changes the energy density by a finite 
amount [5] , however equation predicts a vanishing en¬ 
tropy correction if the operator acts at the boundary of 
the subregion where /3(x) ^ 0. In contrast, backscat- 
tering is known to affect the entanglement entropy most 


strongly at the boundary, and weakly everywhere else [6] . 
This suggests that the entanglement temperature itself 
must be affected by an external field. 

In this manuscript we use AdS/CET duality to show 
that, for a d = 2 CET, a local marginal perturbation 
leads to a local perturbation of the entanglement tem¬ 
perature: P{x) /3{x) + SP{x). Then, from equation 
Q, it follows that 

5S = - j dxdx' l3{x)5P{x'){f{x)f{x')) (3) 

where angle brackets denote expectation value with re¬ 
spect to the unperturbed density matrix. Identifying the 
energy correlation with the specific heat, cy, this equa¬ 
tion is analogous to the Clausius relation of thermody¬ 
namics: 

SS = cvE ( 4 ) 

When /3{x) is associated with projected area of the min¬ 
imal surface on the boundary of AdS, equation @ will 
be shown to express a local relationship between the per¬ 
turbed entanglement temperature and the perturbation 
to the AdS background metric. 

We begin by describing the relationship between entan¬ 
glement temperature and minimal surface area that ap¬ 
pears in the AdS^^+i construction for holographic entan¬ 
glement entropy. Given a minimal surface whose bound¬ 
ary defines a region Q in the CET is, the holographic en¬ 
tropy [7] may be computed from an integral over bound¬ 
ary coordinates: 

The equation of the minimal surface spanning a radius /, 
d—1 dimensional ball of the CET is = z^-fx^Xa, where 
(a = I,..., d—I) are the spatial boundary coordinates. 
Parameterizing the surface {x^ = z{x);x^ = x"}, the 
induced metric is: 

9ab = G%daX^di,X- = + Ab) (6) 
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where = r]y,yl(^ j {fi = 0, — 1) is the un¬ 

perturbed spatial metric for a radius R AdS^^+i space. 
The determinant is ^ = R‘^l^lz^ and, defining the radial 
boundary coordinate = x^Xa^ 

I 27tR 

^ ~ J{^ 

where /3(r) is the inverse of the entanglement tempera¬ 
ture of a d-dimensional field theory, as defined by Klich 
[4] for a radius I subregion of a CFT on an infinite do¬ 
main. Equation 0 expresses that entanglement tem¬ 
perature T^{r) = /3“^(r) at a point r on the bound¬ 
ary is associated with the projected differential area ele¬ 
ment of the minimal surface at the point associated with 
r. Using equation and c = ?>Rl2G^ the holographic 
entropy formula 0 is seen to be completely equivalent 
to integrating the entanglement temperature to find the 
entropy—equation 1.9 of reference [4], written in integral 
form here for d = 2, 


TTC f dx 

3 in/3(a;) 


(8) 


We now consider the effect upon the entanglement tem¬ 
perature, l//3(x), of applying an external field 0o(^) cou¬ 
pled to an operator 0{x) at a point x in the boundary 
theory. Specializing to d = 2 the action of the boundary 
theory is deformed by: 

5" = y* d^x(j){){x)0{x) (9) 

Following the AdS/CFT prescription, the applied exter¬ 
nal field, 00is a source for a scalar field 0(^, x) in the 
gravity bulk, determined by the bulk-boundary propaga¬ 
tor for AdSs {B is a normalization factor): 

f z"^ 

( 10 ) 

where A is the scaling dimension of 0{x) and 5 is a 
normalization constant. We have also written space/time 
coordinates of the boundary theory explicitly. 

(j){z^ x) is a source for the Einstein equations governing 
perturbations to the AdS metric. These perturbations, 
in turn, induce local perturbations to the minimal sur¬ 
face area, thus perturbing 0(x) through equation 0 . As 
seen in the solution ( [1Q| ), a point source at the bound¬ 
ary diffracts as the field propagates in the bulk, effecting 
a spatially extended perturbation of the entanglement 
temperature. This scheme is depicted in figure 0 . 

To confirm this heuristic analysis we have developed 
a procedure to numerically compute the spatial depen¬ 
dence of the entanglement temperature 1/0(x) in a lat¬ 
tice model of one dimensional free spinless fermions de¬ 
formed by a marginal perturbation. In the continuum 
limit such a model realizes a two dimensional c = 1 CFT. 



FIG. 1. Depiction of defect and holographic entropy con¬ 
struction of minimal surface. Localized, marginal scalar per¬ 
turbation of CFT (defect marked by x) propagates into bulk 
AdS (shaded region.) Bulk scalar field sources perturbation 
of metric over an extended region of minimal surface (dashes). 
Perturbation to the minimal surface area projects back onto 
boundary, perturbing entanglement temperature over an ex¬ 
tended region (dashes). 



FIG. 2. Numerical calculation of lattice inverse entanglement 
temperature, 0^, of an 21 = 80 site subregion on a 450 site 
lattice with a defect t' = t/2 at site 57. Unperturbed 0° 
(inset) is indistinguishable from analytic form /3o{x) = 7r(/^ — 
x^)// (not shown) on scale of figure. S/Sx = 0cc — 0S shows the 
singularity at the defect position and algebraic decay (S/Sx ~ 
1/x) corresponding to diffraction of a point source in AdS 
space. 


Our basic numerical results are summarized in figure 
([^. 0^, computed without the perturbation, is nominally 
a parabola, corresponding closely to the discrete version 
of 00(:r) = 7r(/^ — x^)//; the perturbation, ^0^,, has odd 
parity about the defect position and decays algebraically 
corresponding to the spatially extended scheme described 
above. 

We now turn to the calculation of entanglement en¬ 
tropy from SP{x) to arrive at our central result, equation 
Following that, we will work out the relationship 
between 6/3{x) and the perturbation in AdS bulk metric 
that follows from the Ryu-Takaynagi formula for holo¬ 
graphic entropy and gravitational perturbation theory. 
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Writing the perturbed reduced density matrix 

A _ ^Lp- I -1 (/3o(x)-\-dl3(x))f(x)dx 


Equation (18) will be used to relate 5P{x) to the per¬ 


( 11 ) 


turbed AdSs bulk metric. 

The comparable equation to (18) on the gravity side 


with Z set to normalize trp = 1, we would like to ex¬ 
pand the entropy S = —trplogp to first order in S/3{x). 
Let A = f_j^dxP{){x)T{x) and B = J^^dxSP{x)T{x) and 
consider an expansion in /3o{x) keeping only 0{6/3{x)) in 
the exponent: 




gel 


The operator T{x) satisfies the Virasoro algebra 


.TT 


[f{x),T{x')] = D{x, x')f{x) + i^cd^S{x — x') 


( 12 ) 


(13) 


where D{x^x') = 2'Ki5{x — x')dx — 4:7ridxS{x — x') and 
the nested commutators are all linear in T{x). Using 
equation (12), 


relates the perturbation in holographic entropy to the 
perturbed bulk metric, computed from the back reaction 
of the bulk metric to the GET perturbation (|9). At this 
point we also restrict our analysis to time-independent 
marginal perturbations (A = 1). The metric incorporat¬ 
ing a perturbation x) (/i = 0,..., d) of pure AdS 

in the Eefferman-Graham gauge is: 

d 2 

ds^ = -^{dz^ + + h^i,)dx^dx'') (19) 

The resulting entanglement entropy may be found by per¬ 
turbing equation ^ using the specific parameterization 
for the AdSs minimal surface spanning an interval [—/,/] 
on the boundary: x{a) = / sin(cr), 2 :(<j) = / cos(<j) for 
a G [—7r/2,7r/2]. Writing Sg in terms of hxx^ 


P = 


1 - r 




( 20 ) 


J fdo{Xl)SP{X2)D{xiX2)f{xi) 

J ^o{xi)(3o{x2)S/3{x3)D{x2X3)D{xiX2)f{xi) + ... I 

where constant terms arising from the conformal anomaly 
(the last term in (p!^) have been suppressed. Adopting 
the notation 

(...)= tripo ...) = Ttr(e- (^ 5 ) 

where Zq = tr exp (— J Pq{x)T{x))^ the expectation value 
of any operator {0{x)) supported in the interval [—/,/] 
yields its zero physical temperature, infinite domain ex¬ 
pectation value. Making the same expansion of Z and 
noting that {T{x)) = 0, Z may be replaced by Zq in 
equation (HI- Equation may now be viewed as a 
high temperature expansion (in /3o{x)) for each order of 
6l3{x). Similarly expanding the logp factor appearing in 
the entropy definition, 

S = logZo — J dxdx'(3o{x)S(3{x'){T{x)T{x')) (16) 

to lowest order in Po{x) and SP{x), establishing eq. §. 

Eor a GET expressing both chiral modes, the energy 
correlation function regularized by a spatial cut-off a is: 


where we have made use of the relation between central 
charge, c, and properties of AdS: c = 3R/2G. 

Gomparing equations (18) and ( [^ , both expressing 
the central charge, c, suggests a relationship between 
perturbed entanglement temperature and the perturbed 
bulk metric. Gonsider a defect located at boundary co¬ 
ordinate X = 0 and an asymmetric subregion ^ 

centered at x = ^. The perturbed entanglement tem¬ 
perature now depends upon ^ for fixed 1: S/3{x,^). The 
minimal surface may be regarded as a probe in that the 
metric perturbation hxx{x,z) at a point {x,z) may be 
constructed from a numerical computation oi 5p{x,^{z)), 
where ^{z) = x — y/P — z‘^ describes a particular minimal 
surface passing through {x,z). Then, 


hxx{x,z) = ^^^^S/ 3 {x,^{z)) 


( 21 ) 


To compare our numerical computation of the pertur¬ 
bation of the lattice entanglement temperature, 6f3x^ to 
gravitational perturbation theory, hxx is computed for a 
localized applied field in the boundary theory (l)o{x', uj) = 
(j)Q5{x')5{uj). The scalar field in the bulk is found from 
equation ([1Q|) to be: 


(p{z,X,t) = (po 


s/Tk \fAA-~A 


( 22 ) 


This field now sources the Einstein equations through 
the reduced stress-energy tensor: 


{nx)f{x')) = -A 


+ 


2 (x — x' -h iaY {x — x' — iaY 


Performing the x integration in ([l^, 

-V 


AttP 

5S = -^ I dx^ - 


{P — 


(17) 


(18) 


Ttt — + {dx(t>)‘^] 


(23) 


Eollowing reference [3], hxx = hu^ and the linearized 
Einstein equation, 

^zhxx ^^z^XX — 


(24) 
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may be integrated to yield the perturbed AdSs metric, 

= ( 25 ) 

where the external field, 0o and other constants have 
been absorbed into B'. 

In comparing ^{z)) and hxx{x^ z), we focus on the 
important difference in the analytic structure describing 
the defect. Gravitational perturbation theory predicts a 
nonsingular behavior linix^p = B'/4 whereas 

our numerical calculations indicate a singular behavior: 


Sj3(x) ^ —a- 




X — p 


-sign(p) 


(26) 


Numerically integrating equation (|2Q|) using (25) to de¬ 


termine the entropy leads to a nearly constant entropy as 
the defect is moved to the boundary (figure |^. This be¬ 
havior disagrees qualitatively with the exact computation 
of entropy shown in figure ^ which has a pronounced 
trend of large (negative) correction to the entropy as the 
defect approaches the boundary. In contrast, the alge¬ 
braic singularity (26) when integrated in equation (18) 


yields 6S ~ a log i-^, agreeing qualitatively with the ex- 
act entropy calculation as shown in figure ^ (see note 
0 )- The odd parity of SP —contrasting the even parity 
of hxx —appears to be the correct qualitative behavior 
as well: a strong defect that isolates one smaller region 
of size s has an entanglement temperature, Te, that di¬ 
verges as s ^ 0. Therefore, SP < 0 on the side of the 
smaller region and (5/3 > 0 on the side of the larger region. 

The disagreement between gravitational perturbation 
theory and entanglement calculations suggests that grav¬ 
itational perturbation theory may miss a coordinate sin¬ 
gularity associated with the exact dual of this deformed 
CFT. The singular behavior suggested by our calculation 
may bear some relation to the AdS/Schwarzschild solu¬ 
tion used in [8] to study entanglement dynamics following 
a quench. 

We briefly describe our numerical calculation. Us¬ 
ing slightly unconventional notation, our tight-binding 
Hamiltonian is 


= + (27) 

where the kinetic energy density operator is written in 
terms of the single particle kinetic energy matrix and 
fermion operators satisfying the usual commutation rela¬ 
tions. Tx^y = Kx,y^x^y ^x,y — t(,^x,y—l 4“ ^x,y-\-l^' 

(Please note that Roman subscripts now refer to lattice 
coordinates of the tight-binding model.) The marginal 
perturbation, for non interacting fermions, is a weak 
back-scattering term of strength t' acting at a single link 
located at point p. 

The lattice entanglement temperature l/^x is deter¬ 
mined by computing the reduced density matrix for an 
[—/,/] sized sub-region, which Kilch [4] has shown to be 



FIG. 3. Perturbation to entropy as a function of distance 
of defect from boundary (/ — p). Symbols: exact numeri¬ 
cal calculation on a 350-site lattice with a subregion of size 
21 = 160, t' = t/2. Dashed curve: gravitational perturba¬ 
tion theory, equations ( [20| and (25) with B' = 1.6//. Solid 
curv e: boundary theory calculation using equations ( |18| and 
(26) with a = 0.13; see [9]- 


of the form. 


p = Z 5i;/3xTx,y ^28) 

Using the general scheme developed by Peschel uni, it 
may be shown that 


= [log(G-i - /)],,, 


(29) 


where Gx,y = (clcy) is the one particle correlation func- 
tion calculated in the sub-region x,y E [—/,/]• Elements 
of the right hand side of equation ( [2^ that he off the 
tri-diagonal structure of Kx^y were checked for exponen¬ 
tial decay away from the diagonal and discarded. Since 
the eigenvalues of G are exponentially close to 0 or 1, 
sufficient precision must be kept to determine Px- Holo¬ 
graphically, all eigenvalues are needed to resolve the low¬ 
est energy features deep in the AdS space, where z ^ 1. 

However, some further justification is required to inter¬ 
pret these results. In general for small central charge (in 
our case c = 1) or small N SU(A^) GFTs, quantum cor¬ 
rections to the Ryu-Takayanagi (RT) formula are needed 
in the bulk to compute quantities other than the entan¬ 
glement entropy. For instance, according to the semi- 
classical RT formula, the mutual information of well sep¬ 
arated regions approaches zero, and quantum corrections 
give the dominant non-vanishing contributions m- Sim¬ 
ilarly, computation of the Renyi entropies in the gravity 
dual do not have a simple semiclassical geometric inter¬ 
pretation. Since the density matrix captures the entire 
operator content of the theory, it is natural to question 
whether a semiclassical gravity computation (or pertur¬ 
bation thereof) can fully describe the density matrix of 
the boundary theory [12]. 

The c = 1 GFT used in our calculations can be realized 
as free fermions for (Luttinger parameter) ^ = 1 or com- 
pactified free bosons for arbitrary g. In a noninteract¬ 
ing theory, the reduced density matrix takes the simple 
quadratic second quantized form in equation (28) and the 
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semiclassical geodesic appearing in the RT construction 
of entropy simply encodes the entanglement temperature 
function, ^{x). Specifically, the equation of the geodesic 
is 


Z^{x) = —ip{x) 

TT 


(30) 


In this special case, the semiclassical physics of the 
bulk completely describes the boundary theory density 
matrix for a single subregion. When the boundary theory 
is deformed by a marginal perturbation, we have shown 
numerically that the quadratic form of the density matrix 
is preserved, but with a perturbed entanglement temper¬ 
ature and entropy. On the gravity side, the marginal 
perturbation sources the AdS space in the neighborhood 
of the RT geodesic, perturbing the entropy. It is reason¬ 
able to expect these two computations to be comparable, 
and the fact that the central charge appears identically in 
both computations supports this conclusion. In this way, 
for the particular case of c = 1, we expect perturbations 
to the bulk AdS metric in the vicinity of the minimal sur¬ 
face may be extracted from the boundary computation 
of the entanglement temperature. 

In this manuscript we have broadened the thermody¬ 
namic analogy that presently links quantum information 
and gravitational physics by uncovering a generalization 
of the Clausius relation. It is shown that application of 
an external field, coupled to a marginal operator in the 
boundary field theory, effects a local change in the en¬ 
tanglement temperature. The resulting entropy change 
is then g iven by an equation analogous to the Clausius 
relation with spatial energy correlations at zero 

physical temperature playing the role of heat capacity. 

By comparing this relation to the holographic expres¬ 
sion for entropy, we have identified a local quantity in 
the boundary theory, S/3(x, that appears to corre¬ 


sponds to a local degree of freedom in the dual grav¬ 
ity theory, the perturbed bulk metric hxx{x^z). We 
are therefore able to construct perturbations in the bulk 
geometry from direct diagonalizations of the boundary 
quantum theory. For the particular quantum impurity 
model we study, there is a significant difference between 
the analytic structure of the metric determined computa¬ 
tionally from the boundary theory—it is singular at the 
position of the defect—and the metric calculated by grav¬ 
itational perturbation theory. In future work, it will be 
important to study extended perturbations, as opposed 
to a localized one, in that our computations suggest dra¬ 
matic perturbations to the bulk geometry. 


We note that equation (§ may be generalized to fi¬ 
nite physical temperature by using the appropriate P{x) 
from reference [4] and crosses over to the thermody¬ 
namic Clausius relation at temperatures larger than 
the corresponding inverse length scale (or mass scale 
in the case of a system that is gapped). Lastly we 
add that SP{x) might be calculated analytically from 
tr(pO(x)) = {0{x)) where the left hand side is expanded 
to 1st order in 6^ (as described leading to equation (14)) 
and the right hand side is an ordinary perturbative ex¬ 
pansion to 2nd order in f /t at zero physical temp. This 
is actually the basis of Peschel’s numerical scheme in 
[To]; however, we have only been successful in computing 
S/3{x) numerically. 


We wish to point out that recent work has come to 
our attention that studies gravity duals of marginally 
deformed CFTs from a quantum information perspec¬ 
tive [13]. Our research was supported by Research Cor¬ 
poration CC6535 (CL) and Howard Hughes Medical In¬ 
stitute Scholar Program (BC). GL wishes to acknowl¬ 
edge many useful discussions with Benjamin Burrington, 
Adam Durst, Barry Friedman and Daniel Yee. 
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